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The impact of 1 GeV/Nucleon U^'^^ projectiles against atomic hydrogen is studied by direct 
I numerical resolution of the time-dependent non-relativistic wave-equation for the atomic electron 

^— »' . on a three-dimensional Cartesian lattice. We employ the fully relativistic expressions to describe 

the electromagnetic fields created by the incident ion. The wave-equation for the atom interacting 
, with the projectile is carefully derived from the time-dependent Dirac equation in order to retain 

^ ■ all the relevant terms. 
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Atomic ionization by an electromagnetic wave, either incoherent or coherent, has been widely studied since the old 
times of Quantum Mechanics in many different situations. It is also obvious, but much less common in this context, 
to study ionization by the electromagnetic fields generated by a fast charged projectile passing nearby. In the case of 



^ a laser field, or any other electromagnetic wave, the field is a radiation field far from the sources, in other words, it is 
, made of real photons 0]. However in case of a rapidly passing projectile, the field has different properties and it is 
^ ' just made of virtual photons. 

That those virtual photons can ionize the atom is completely clear, but the ionization dynamics should be quite 
different from the laser field case. Recently, considerable interest has been devoted to this subject ||^, both for single 
and multiple ionization of atoms with fast ions. However, very few experiments have been reported with relativistic 
incident ions. At the GSI in Darmstadt experiments with relativistic incident ions have been performed. The 
' GSI experiment studied the single and double ionization of helium by 1 GeV/Nucleon U^^+ impact. As these authors 
state, the relativistic ion generates a sub-attosecond super-intense electromagnetic pulse. 
I Due to the geometry of the system (incident projectile and nucleus), there are no symmetries present except for a 

I mirror symmetry along the plane defined by the projectile trajectory and the initial position of the nucleus. Therefore 
. three-dimensional studies are needed and this is extremely difficult to be done ab initio for a two electron system. 

In the present paper we present a very realistic description of the (relativistic) ion-atom interaction but just for one- 
, electron atoms. We compute the electron wave-function in a Cartesian three-dimensional lattice, taking into account 
' the relativistic dynamics of the projectile. 
^ , Calculations in three-dimensional Cartesian grids have been shown to be a very good description for the p + H 
' collisions with non-relativistic incident protons. The agreement with experimental results is remarkably good 
^ _ for the range of parameters considered. Collisions with antiprotons have been also investigated with these ab initio 
simulations However these calculations are done in a context that is not able to describe the peculiar features 

'■q^' of the electromagnetic fields (electric and magnetic fields) generated by the relativistically moving projectile. To our 
• • > knowledge no previous studies have been done based on the numerical direct resolution of the wave-equation, in cases 
^ ' where the incident ion is relativistic. 

, A remarkable side consequence has been found. We need to describe a Lorentz transformed Coulomb field instead 
5^ ' of the typical radiation fields studied in photo-ionization. For the correct quantum description of such interaction, it 
^ [ is necessary to include one extra term in the time-dependent Schrodinger equation. This term is missing in most of 
the studies on laser photo-ionization, that start from a time-dependent Schrodinger equation conceptually wrong. 

The paper is organized as follows. In Sec. II we describe the electromagnetic fields created by a relativistically mov- 
ing ion. In our calculations, we neglect both the recoil of the hydrogen nucleus so as the change in the trajectory of the 
ion: these approximations are justified in terms of classical simulations (Monte Carlo) in Sec. HI. The derivation and 
discussion of the time-dependent wave-equation for the atomic system is described in Sec. IV. A detailed description 
of the integration algorithms and numerical techniques is also given. In Sec. V we establish the parameters of the 
calculations and we present the obtained numerical results. The final Section is devoted to conclusions. 
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II. ELECTROMAGNETIC FIELDS 



As our particular choice of the coordinates, to simphfy the notation without loss of generality, we consider that 
the heavy ion projectile is moving in the plane z = in parallel direction to the x-axis. The hydrogen nucleus is 
initially at the origin and the projectile trajectory is Xion{t) = xq + vt = Xo + Pet, yion(t) = b and Zion{t) =0. 6 is the 
impact parameter (minimum distance to the nucleus in the trajectory of the ion). The notation P = v/c is standard 
in relativity. The projectile is so heavy (Mio„ — 238mAr, being mjv the nucleon mass) and so energetic that we can 
perfectly assume that its trajectory is not affected by the atom. This point is discussed in Sec. Ill in terms of classical 
simulations. 

It is not difficult to show that the electric and magnetic fields created by such a relativistic projectile at an 
arbitrary point (x, y, z) of space are. 



E(r,t) = - 



Z 67 



and 



B(r,i) 



2 

7^ (a; — xo — vt) + {y — 6)^ + 

— Ze7^ v/c 
7^ {x — xq — vt)^ + {y — b)'^ + 



\{x -xq- vt) , y, z] 



(1) 



3/2 



(0,z,-y) 



(2) 



with 7 = l/-y/l — Z indicates the projectile ion charge, that for the case studied here is Z=92. 
These fields are generated by a scalar potential 



and a vector potential 



77(r,i) 



A(r,i) 
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\J j'^ (x — xq — vt)^ + {y ^ b)^ + z^ 
Z ejv/c 



\J ^"^ {x — xq — vt)^ + (y — bi)^ + z"^ 
The velocity vector has only one component, v = (w, 0, 0). Therefore, A(r, t) = Ax (r, t) Gx and Ay = Az = 0. 



(3) 



(4) 



III. CLASSICAL SIMULATION 



Although we give, in this paper, a quantum description of the dynamics of the atomic electron interacting with the 
ion, it is worth to start with a classical simulation of the three-body problem studied. For the considered parameters 
of the incident ion, it is almost trivial to perform classical Monte Carlo simulations for the motion of the projectile 
ion, the proton and the electron. The goal of such simulations is just to introduce the right approximations for the 
subsequent quantum description. 

A few conclusions appear from the classical simulations for heavy projectiles (U^'^+) that move relativistically (1 
GeV/Nucleon, or more). First, the projectile trajectory is not changed at the space scale and with the precision 
that we are interested in. Second, the hydrogen nucleus is accelerated after the collision with the projectile but the 
momentum transfer is negligible unless we consider impact parameters smaller than 0.1 a.u. We are not interested in 
head-on collisions, that could be relevant for nuclear physics. 

In Fig. we represent the motion of the hydrogen nucleus for different values of the impact parameter (6). The 
horizontal axis shows the motion in the longitudinal direction (x coordinate), and the vertical axis shows the motion 
in the transversal direction (y coordinate), where the impact parameter has been added for better observing all the 
curves. Both scales are in atomic units. The trajectories has been plotted since the incident projectile is placed at 
x{t) = — 10 a.u. (marked with big circles) until it reaches x{t) = 10 a.u. 

A result is clear: for the study of the ionization up to the level of accuracy we are interested in, it is fairly exact to 
consider that the projectile ion follows its trajectory unaltered. It is also reasonable, due to the short interaction time, 
to consider that the hydrogen nucleus remains unaltered at its initial position. Therefore only the electron motion 
needs to be described in detail. On the other side, the classical simulations indicate that the electron's motion is 
essentially non-relativistic, except for extremely small impact parameters. 
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IV. TIME-DEPENDENT WAVE-EQUATION 



Because the electron is interacting with such time and space dependent scalar and vector potentials, it is worth 
to carefully describe the structure of the wave-equation. It is well known that the Dirac equation for an electron 
interacting with an arbitrary electromagnetic field E(r, t) and B(r, t) described with the scalar potential, ry (r, t), and 
the vector potential, A (r,t), can be written as a second-oder differential equation |^ 

d \^ 

ih— + er] + eVj $_d = (cp + eA)^ $d + rn^c'^ + ehc cr' • B <I>d - iehc a ■ E , (5) 

where is the standard Dirac four-spinor, and cr' and a are 4x4 matrices. The space and time dependence on 
the fields, potentials and wave-function has been removed for simplicity in the expressions. To describe the Coulomb 
potential of the hydrogen atom we have also included the eV term. 

To move to the non-relativistic domain we start by neglecting the —iehc a • E $£) term, as it mixes particle and 
antiparticle states. We also neglect the coupling between spin and magnetic field, cr' -B, because the degree of freedom 
of spin is not expected to play an important role in the dynamics of the electron, in the cases that we are considering. 
Now we can just consider a one-component wave-function $ 

i^TT +eri + eV] $ = (cp + eAf $ + m^c'' $ . (6) 
ot J 

Expanding the square that includes the time derivative, and taking into account that the scalar potential can be 
time- dependent: 

-h^^+ 2iehir] + V)^ + leh ( ^) $ + e^irj + Vf<S> = (cp + cA)' $ + m^c^ $ . (7) 

o''t ot \<yt J 

We now introduce explicitely the fast time oscillation, 

$(r,t) = 1'(r,t) expf-z— ij (8) 



and the wave-equation becomes. 



9^ 

2ihmc^— - fi^-g^ + 2emc^{rj + V) + 2i h eirj + V)-g^ 

+ieh^ + e^{r] + Vf = {cp + eAf ^ (9) 

so that we have eliminated the fast oscillation due to the mass term. The next step towards the non-relativistic 
wave-equation consists of neglecting the d^'i'/d^t term, because the fast oscillations due to the mass term have 
been explicitly accounted for. Moreover it is consistent with the non-relativistic limit we are looking for to assume, 
mc^ >> e{ri + V). Thus the wave-equation is simplified to, 



2ihmc — — 
dt 



2emc^(ri + V) + ieh—^ 
ot 



* = (cp + cA)^ (10) 



The iehdr]/dt term should not he eliminated because even when the modulus is very small compared with the electron's 
rest mass, it is dephased due to the i factor. Therefore it is necessary to keep this term, 

^ d'if 1 / e \ 2 eh dri , , 

dt 2m V c J ' 2mc2 dt ^ ' 

It is worth to write this wave-equation in the form of a Hamiltonian system. This Hamiltonian has an extra term 
idri/dt, that needs to be included in the time-dependent Schrodinger equation. This term is, however, fundamental for 
the validity of this non-relativistic equation for any scalar and vector potentials with arbitrary space-time dependences. 
Expanding the (p + e/cA)^ term, and taking care of the non-commutability of p and A, one gets, 

^ p1 ^ + ^A • p * + ^^A2 -i^-eri^-eV^f- i— (v-A+ * . (12) 

ot 2m mc 2mc 2mc \ c dt J 
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If we introduce the Lorentz condition for the scalar and vector potentials, 



c ot 



the wave-equation is heavily simplified, 



ih 



dt 



2m 



mc 



-A • p vj/ 



^ - er^ ^ - e-qV . 



(13) 



(14) 



Observe that this equation is different from the standard time-dependent Schrodinger equation, that considers 
the term (cp -I- eA)^: it includes the divergence of the vector potential but it does not include the very important 
idrj/dt term. Such equation is not general and can not be safely employed to describe the interaction with arbitrary 
electromagnetic fields. Only transversal electromagnetic fields can be accurately described. We have introduced this 
derivation prior to start computing the dynamics of the electron in order to clarify this point, that could become a 
source of error in the description of the system. 

One final remark on the above equation. The fields that we describe with A and 77, Eq. (^ and Eq. (^, are invariant 
under Lorentz transformations. The time-dependent Schrodinger equation is invariant under Galilean transformations 
due to its non-relativistic nature. However, Eq. (14) shows neither one nor the other kind of invariance. This formal 
inconsistency in our model is not important if the dynamics of the electron is non-relativistic (like in the cases that we 
are considering), and thus Lorentz transformations reduces to Galilean transformations (in the low velocity regime). 

The final equation that we consider is: 



ih 



dt 



V^*(r, t) - i — A(r, t) • V *(r, t) 

2m ^ ' mc ^ ' ' 



2mc' 



+ 



:A(r, tf *(r, t) - eri{r, t)*(r, t) - eV{r) ^'(r, t) . 



(15) 



V{r) — e/ x'^ + y'^ + z'^ is the Coulomb potential. The electromagnetic potentials are given by Eqs. and (Q), and 
they already verify the Lorentz condition because they come from the Lorentz transform of a Coulomb potential. 

We numerically solve Eq. ( p^ on a three-dimensional Cartesian lattice with a uniform grid spacing Ax = Ay = 
Az — 0.3 a.u. The lattice is taken in such a way that avoids the singularity in the origin by considering grid points 
at half integer times the grid spacing. 

Our numerical technique is based on a symmetric splitting of the time-evolution operator so that the wave-function 
at a given time t + At is calculated from the wave-function in the previous time step t: 



5'(r,t At) ~ exp 
~ exp 
X exp 
X exp 



'i—H(t + At/2) 
h 

'i^H,{t + At/2) 
'i^H4t + At/2) 
-^^H^{t + At/2) 



*(r,i) 
exp 
exp 

L z,ii 
^'(r,i) + 0(At3) 



.At 

.At 
W 



^7^Hy{t + At/2) 



^^Hy{t + At/2) 



(16) 



We have defined the hamiltonian operators as follows 

f,2 Q2 „2fe2 



d 



= -TT-T^ + 7r-,Mr,ty - ^_^,(r,^)— - e^{r,t) - eViv) 



2m dx^ 
2m dy"^ 



2mc 



mc 



dx 



2m dz 



(17) 



The exponentials are thus expressed in the Cayley unitary form 0, that preserves the norm of the wave-function, 
and a Crank-Nicholson scheme is employed to evaluate the space derivatives. Similar methods have been successfully 
employed to integrate the non dipole Schrodinger equation for the interaction of atoms with very intense laser fields 



The initial state for our calculations (the ground state (/'is(r) of atomic hydrogen) is computed in a 100 x 100 x 100 
lattice with the imaginary time propagation method. The computed energy for this state is Eb = —0.490 a.u. 
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Two different grid sizes have been employed to describe the ion-hydrogen coUision, depending on the value of the 
impact parameter (b): 



N^xNyxN,= 



280 X 380 X 380 for 6 < 10 a.u. 
280 X 280 X 280 for 6 > 10 a.u. 



(18) 



Absorbers (a mask function) were employed at the integration boundaries in order to avoid reflections of the ionized 
population. The mask function has the form sin^/^ and it is applied over 40 points along the edge of the grid. However, 
the interaction time is short enough to avoid a very large amount of population reaching the boundaries (absorbed 
population is always smaller that 1 percent). 



In our calculations, the hydrogen nucleus is placed at the origin of the Cartesian coordinates. For the selected 
energy of the incident ion (1 GeV/Nucleon) the relativistic parameters are p = 0.3696 and 7 = 1.076. The U^^"*" 
nucleus is initially at 2;(0) = a;o = —200 a.u., y(0) = b and z(0) = a.u. That choice of the initial value of a;(0) is 
a good compromise between accuracy -due to the long range of the Coulomb potential- and a reasonable computer 
time. The electronic wave-function is propagated in time until the ion reaches x{t) = 600 a.u. (t = 15.8 a.u. is the 
final time). 2667 time iterations are employed, so that the ion covers a distance equal to the grid spacing (0.3 a.u.) 
each time step {At — Ax/v — 0.00592 a.u.). In order to avoid problems with the singularity of the Coulomb potential 
of the ion, the projectile moves equidistantly to the nearest points of the lattice. 

We investigate the effect of the impact parameter in the dynamics of the atomic electron by varying b (non- 
uniformly) from 1 a.u. to 150 a.u. We avoid smaller impact parameters because they contribute negligibly to electron 
ionization cross sections and they turn out to be relevant for ion-proton scattering (a physical situation out of the 
scope of the present paper). 

We evaluate, during the time propagation, the projection of the wave function over the ground state '/>is(r) to 
compute the total excitation probability (probability of finding population both in free states as in bound excited 
states), for each value of the impact parameter: 



V is the volume of the integration box. 

In Fig. ^ it is represented the computed values of Pexc for different impact parameters b, at the final time t = t. 
In fact, due to the fast motion of the ion, very few iterations after the ion has passed at the minimum distance to the 
nucleus {b) , the excitation probability remains approximately unchanged. From these values it is possible to estimate 
the total excitation cross section at the final time, defined as 



We interpolate the excitation probability at the intermediate points with cubic-spline standard methods. The obtained 
value is aexc{T) ^ 8.4 x 10~^^ cm^. 

To gain insight into the dynamics of the electrons, we compute the population that remains, at any time of the 
interaction, at a distance R to the nucleus greater than 15 a.u. [PR>i5(i)] or greater than 20 a.u. [PR>2o{t)], that 
can be regarded as ionized electrons. At the final time t, both values of the population have been plotted in Fig. ^ 
in terms of the impact parameter b. Dashed lines and diamonds corresponds to Pr>i^{t) while solid lines and circles 
correspond to PR>2a{T). Both curves sharply decreases from 6 = 1 to & = 5 a.u., so that very few electrons can be 
measured for larger impact parameters. 

In Fig. ^ we plot the time evolution of the population PR>i5{t), for different values of the impact parameter. The 
horizontal axis contain the interaction time in atomic units and the vertical axis the population PR>i5{t). We may 
approximately know the velocity distribution of the ejected electrons, by simply assuming that ionization process 
occurs at the time when the ion is placed at the minimum distance to the nucleus [xioniU) = 0, yion{ti) = b and 
ZioniU) =0]. In fact this is a good approximation due to the very high velocity of the projectile. The time ti at 
which the ion is placed at this minimum distance is marked with a vertical arrow. Vertical dashed lines indicate the 
minimum necessary velocity of the ejected electrons, Vm, to reach the distance R — 15: Vmit) ~ (t — ti)/15 a.u. (for 
t > ti). Observe that no electrons are ejected with velocities greater that 5 a.u. {v = 0.037c). Thus, to employ the 
non-relativistic wave-equation is justified to describe the electron dynamics. 



V. NUMERICAL RESULTS 




(19) 




(20) 
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At the final time r, we have also computed the angular distribution of the electrons. To do so we consider electrons 
placed between two spheres, i.e. electrons placed at a distance R from the nucleus such that Ri < R < R2 with Ri = 15 
a.u. and R2 — 30a. u. These values of the radii mean that we are considering electrons ejected with approximated 
velocities in the range 1.3 < Vm < 2.6 a.u. In Fig. |^ this angular distribution has been plotted for the impact 
parameter b = 1 a.u. (plots at the top) and 6 = 3 a.u. (plots at the bottom). The angular distribution per volume 
unit is defined by: 

(21) 



where we have defined: 



n{Ri < R< R2,0,< 



pR2 pe+A0/2 ^0+A0/2 

P{Ri < R< R2,e,<l)) ^ / / \^{r,e,(j))\^r^ sin 0drd0d(f> 

I-R2 i-e+Ae/2 ^0+A0/2 

^{Ri<R<R2,0A)^ I I I smOdrdedcj) . (22) 

J Ri Je-A0/2 J0-A0/2 

The spherical coordinates are defined in the standard form, x = r sin cos t^, y = r sui0 sin cj) and z = rcos0. The 
angular step is A0 = A(/) = ir/lS rad. Plots on the left are contour plots with dark regions representing the maximum 
population. Contour lines are in linear scale. Plots on the right are surface plots, also in linear scale. 

Observe the clear difference between both pairs of plots. For the impact parameter 6=1 a.u., the peak of the 
ejected electrons is placed along the positive part of the y-axis {0 — tt/2 and (p — 7''/2) but also appears a second 
peak along the negative part of the y-axis (9 = tt/2 and </> = 3tt/4). This second peak is a consequence of the fact 
that the ion crosses well over the electronic cloud. The minimum density of electrons are found along the x-axis, the 
longitudinal direction, {9 — 7r/2, = and = 7r/2, (j) = tt). However for the impact parameter 6 = 3 a.u. only the 
first peak can be observed: in this case the ion does not cross over the electronic cloud. 

Finally, in Fig. |^ we plot the probability density at different times of the interaction with the ion, at the planes 
z = (plots at the top), y = (plots in the middle) and a; = (plots at the bottom). The impact parameter is 
6 = 3 a.u. The column on the left corresponds to the time at which the ion is placed at the minimum distance to 
the nucleus [xionit) = 0]. The column in the middle is for Xionit) = 400 a.u. {t = 11.8 a.u.) and the column on the 
right for Xionit) = 600 a.u. (t = t = 15.8 a.u.). The contour lines are in logarithmic scale. In the frame z = in the 
first column, the effect of the incident ion can be observed as a distortion of the wave-packet. Once the projectile has 
passed, ionized population moves mainly along the transverse direction [y axis). The slight trembling in the lower 
contour lines that can be appreciated in a few frames, is due to spatial discretization effects and the non continuous 
evaluation of the trajectory of the ion. 

The same arrangement of the plots is used in Fig. but now for the impact parameter 6 = 10 a.u. The distortion 
of the wave packet is similar to that shown in Fig. However, ejected electrons are even slower and at the final time 
the spreading of the wave-function is not so marked. 



VI. CONCLUSIONS 



We have presented very realistic three-dimensional numerical simulations to describe the collision of relativistic 
heavy ions against hydrogenic atoms. Our simulations are based on the ab initio integration of the time-dependent 
wave-equation derived from the non-relativistic limit of the Dirac equation, that is valid for arbitrary electrogmanetic 
potentials. This equation differs from the standard time-dependent Schrodinger equation in a term that includes 
the time derivative of the scalar potential: neglecting this term leads to a conceptually wrong formulation. The 
recoil of the ion and the motion of the atomic nucleus have been neglected in our simulations. We have employed a 
three-dimensional Cartesian grid to integrate the wave-equation. As results, we should point out that the probability 
of excitation of the atom sharply increases for impact parameters 6 < 15 a.u. However, the ejection of fast electrons 
is only relevant for even smaller values of the impact parameter (6 < 5 a.u.). The space distribution of such electrons 
is centered along the transversal direction. It is composed of one maximum (in the positive direction of the y axis) 
even for short impact parameters (6 > 3 a.u.). Only in collisions where the ion crosses over the electronic cloud a 
more complex structure appears. 
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FIG. 1. Classical trajectories of the atomic nucleus for different values of the impact parameter b. The horizontal axis 
represents the motion in the x direction and the vertical axis represents the motion in the y direction, where the impact 
parameter has been added for better observing the different trajectories. Distances are in atomic units. 
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FIG. 2. Excitation probability Pexc [see Eq. (|l9|)] at the final time t — r, in terms of the impact parameter b. 
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FIG. 3. Probability of finding population at distances greater than 7? = 15 a.u. from the nucleus (Pr>ib, in dashed lines 
and diamonds) and greater that R = 20 a.u. (Pji>20 in solid lines and circles), at the final time t = t. 




Time (atomic units) 



FIG. 4. Population at distances greater than _R = 15 a.u. during the interaction time (Pii>i5(t)), for impact parameters 
from 1.2 to 10 a.u. An approximated spectrum of the ejected electrons can be estimated for different values of the impact 
parameter. See text for explanation. 
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FIG. 5. Angular distribution of the ejected electrons for the impact parameter 6=1 a.u. (plots at the top) and 6 = 3 a.u. 
(plots at the bottom). The plots show the population that at the final time is contained between 7i = 15 a.u. and R = 30 a.u. 
Contour lines are in linear scale, with dark regions indicating more populated regions. 
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FIG. 6. Probability density at different times of the interaction with the ion, at the planes z = (plots at the top), y = 
(plots in the middle) and x = (plots at the bottom). The impact parameter is 6 = 3 a.u. The column on the left corresponds 
to the time at which the ion is placed at the minimum distance to the nucleus [xion{t) = 0]. The column in the middle is for 
Xionit) = 400 a.u. (t = 11.8 a.u.) and the column on the right for Xion{t) = 600 a.u. (t = t = 15.8 a.u.). The contour lines are 
in logarithmic scale. Atomic units axe used. 
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FIG. 7. Probability density at different times of the interaction with the ion for 6 = 10 a.u. The arrangement of the plots is 
the same as in the previous figure. 
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